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Negative Poisson's ratio and semi-soft elasticity of smectic-C liquid crystal elastomers 
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University of Surrey, Guildford, GU2 7XH, United Kingdom. 

Models of smectic-C liquid crystal elastomers predict that it can display soft elasticity, in which 
the shape of the elastomer changes at no energy cost. The amplitude of the soft mode and the 
accompanying shears are dependent on the orientation of the layer normal and the director with 
respect to the stretch axis. We demonstrate that in some geometries the director is forced to rotate 
perpendicular to the stretch axis, causing lateral expansion of the sample; a negative Poisson's 
ratio. Current models do not include the effect of imperfections that must be present in the physical 
sample. We investigate the effect of a simple model of these imperfections on the soft modes in 
monodomain smectic-C elastomers in a variety of geometries. When stretching parallel to the layer 
normal (with imposed strain) the elastomer has a negative stiffness once the director starts to rotate. 
We show that this is a result of the negative Poisson's ratio in this geometry through a simple scalar 
model. 

PACS numbers: 83.80.Vx, 61.30.Vx, 46.32.+X, 62.20.dj 
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I. INTRODUCTION 

Liquid crystal elastomers (LCEs) are soft solids com- 
posed of flexible polymers, with attached liquid crys- 
talline mesogens, crosslinked into a network pQ. A va- 
riety of liquid crystalline phases of LCEs have been syn- 
thesized, including the nematic and smectic phases. The 
nematic phase undergoes deformation at no energy cost, 
known as soft elasticity [2 in both monodomain sam- 
ples stretched perpendicular to the director [3 , and in 
some types of polydomain samples \^^. Soft elastic- 
ity in nematic elastomers requires several sympathetic 
shears to develop during the deformation as the direc- 
tor rotates. The deformation of the sample must also 
obey the boundary conditions imposed by the clamps. 
Thus, the sample forms a striped microstructure on the 
micrometre length scale consisting of domains in which 
the director rotates in opposite directions in adjacent do- 
mains [T . Without this microstructure soft deformations 
would not be possible. In nematics this microstructure 
has been observed in some detail experimentally [8 , and 
its mathematical properties described [9^. Theoretically 
an ideal nematic LCE should be perfectly soft, however 
in practise a small force must be applied to deform the 
LCE. This semi-soft behaviour is due to various imper- 
fections in the elastomer, and can be incorporated into 
theoretical models by the addition of a semi-soft energy 
term that penalises rotation of the director with respect 
to the rubber matrix [10 . 

Smectic LCEs have been fabricated in both the 
smectic- ^4 (Sm-A) and smectic-C (Sm-C) phases, and 
with both main chain [TT] and side chain [12] architec- 
tures. Their mechanical behaviour can be modelled by 
adding in the embedded smectic layers to the nematic 
elasticity free energy [131 113 • Sm- A elastomers with a 
high degree of smectic order exhibit a sharp change in 
their elastic deformation when deformed parallel to the 
layer normal, and are extremely anisotropic materials, 
behaving as 2-D elastic materials [T2j[T5l[T6]. However, 



the response of smectic elastomers depends on the chem- 
istry, the crosslinking procedure, and domain sizes in the 
sample [17 . In contrast to the nematic phase the Sm-A 
phase does not show any soft elastic behaviour because 
the director is locked parallel to the layer normal. How- 
ever, the elastic behaviour of Sm-C elastomers is pre- 
dicted to be more complex. The director is free to rotate 
on a cone around the layer normal, with fixed tilt angle 
as shown fig. IT] a) . As a consequence it is predicted to 
have a soft elastic mode just as in nematic elastomers 
[TSl [To] . A more complicated combination of shears is 
required in Sm-C soft modes. As a result of the com- 
patibility requirements between these deformations a far 
more restricted set of tensile geometries are predicted to 
deform softly with clamped boundary conditions [20] . 

To test these theoretical results experimentally, a mon- 
odomain must be produced which requires alignment of 
both the layer normal and the director. Using a two-stage 
crosslinking method, the director field can be uniformly 
aligned [21]. However, the layer normals are tilted at a 
fixed angle on a cone around the director. We will refer 
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FIG. 1. a) The layer normal and director in the Sm-C phase 
(polymer chains not shown) , b) the director and layer normals 
in a pseudo-monodomain. 

to this as a pseudo-monodomain (see fig. [lib)) [22] in or- 
der to distinguish it from a monodomain, which has both 
the layer normal and the director uniformly aligned, and 
polydomain, which has random director and layer normal 
alignment. The layer normals in the pseudo-monodomain 



can be aligned by a second uniaxial deformation per- 
pendicular to one of the layer normals ^23^, or alterna- 
tively by a shear deformation perpendicular to the direc- 
tor [24 . To our knowledge no mechanical experiments on 
monodomain Sm-C elastomers have been reported, per- 
haps because of the difficulty in aligning these samples. 
However, the spontaneous deformations associated with 
changing the phase of the elastomer from Sm-A to Sm-C 
have been observed [25 . Many more experiments have 
been carried out on the more accessible polydomain sys- 
tem [TTl [26] . Unfortunately, as with many polydomain 
systems, this is more difficult to model theoretically. 

The director reorientation in soft modes may be par- 
ticularly important for the electromechanical properties 
of chiral Sm-C* phase elastomers [27 . The liquid crys- 
tal rods have a permanent dipole moment that is aligned 
perpendicular to both the director, and the layer nor- 
mal. The coupling between the macroscopic mechanical 
deformations, and the microscopic orientation of these 
dipoles results in their piezoelectric properties. These 
materials show both the direct piezoelectric effect in 
pseudo-monodomains [28], as well as the inverse piezo- 
electric effect [29] . Spontaneous polarization of pseudo- 
monodomains has also been reported [22 . 

This paper is organised as follows. We will describe 
the model of Sm-C elastomers that will be use in ^TTj 
and show that the soft modes in this model have negative 
incremental Poisson's ratio in some geometries. In ^III 
we will illustrate the effect of the semi-soft elastic term 
in four different geometries. We will then summarise the 
effect of this term, and discuss the model predictions in 
relation to the mechanical experiments in polydomains 
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II. MODEL FREE ENERGY 

The model of a Sm-C elastomer that will be used here 
is described in Refs. j30l[3T]. The free energy has contri- 
butions from the nematic elasticity Fnem-, smectic layer 
spacing Fsmi and the energy penalty for changing the tilt 
of the director with respect to the layer normal Ftut- The 
nematic elasticity is given by 



\iiTt A-|o-A^-r^ 



(1) 



where ji is the rubber shear modulus, and A is the defor- 
mation gradient. The step length tensor before the defor- 
mation has been applied is |o = ^+ (r — l)nono, with no 

the initial director, ^ the unit tensor, and r the polymer 

anisotropy. The current step length tensor is denoted 
by |, and its inverse by l""*^ = ^ + (1/r — l)nn, with n 

the final director. In principle a Sm-C elastomer should 
have a biaxial shape tensor for the polymer backbone 
because its shape may be affected by both the director 
alignment and the layer normal direction. For simplicity 
we will approximate it as uniaxial here, depending only 



on the director. We will also assume that the nematic 
and smectic order parameters remain fixed throughout 
the deformation. 

It is assumed that the smectic layers are embedded 
in the rubber matrix, so that the corresponding layer 
normals k will deform like embedded planes 
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where ko is the initial layer normal. The layer spacing is 
penalised by the smectic liquid crystal modulus B 
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where d is the final layer spacing, and do is the initial 
layer spacing, and d/do = 1/|A~"^ • ko|. Fg^ describes 

the free energy penalty for deviations of the layer spacing 
away from that required to accommodate the smectic 
mesogens. For tilted smectic mesogens, the required layer 
spacing is cos^/cos6>o, where the their tilt angle with 
respect to the layer normal is ^o in the initial state, and 
in the current state. The free energy term that penalises 
the deviation of the director from a tilt angle ^o is 



tilt 



^at (cos^ Oo — {n ■ k)^ 



(4) 



where at is the tilt modulus, and n • k = cos 6. 

It will be assumed here that the bulk modulus of 
the rubber is much larger that the shear, tilt and 
smectic moduli, so that the deformation gradient obeys 
det[A] = 1, hence it conserves volume. Typically the 

smectic layer modulus is very large compared to the rub- 
ber shear modulus, i.e. 5 ^ /i (at least in smectic elas- 
tomers of a similar type to that of Nishikawa et al. [15 ), 
so that the layer spacing remains almost fixed. The tilt 
modulus is also large compared to the shear modulus 
at ^ /i, so that the tilt angle remains close to ^o [32j . 



A. Soft elasticity 

The free energy outlined above permits the subset of 
the nematic soft modes that maintain the layer spacing. 
There is only one soft mode that satisfies this (up to a 
global rotation), and it corresponds to a rotation of the 
director about the layer normal [18l |19] . We summarise 
some of the properties of this soft mode here, as they are 
crucial in understanding the semi-soft response of the 
elastomer. 

We assume that the layer normal points along the z 
direction, and the director is tilted into the y direction, 
i.e. ko = z and no = z cos 6>o + y sin 6>o in the starting 
state. The soft modes can be parameterised by the angle 
(/) which gives the rotation of the director no around the 



layer normal towards the x direction. The deformation 
matrix is given by [18 
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The deformation components as a function of rotation 
angle are illustrated in Fig. [2] 
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FIG. 2. For the parameter values r — 2 and ^o = 0.5 radians 

a) shows the diagonal components of the deformation matrix, 

b) shows the shear components, and c) shows an illustration 
of the deformations on the LCE, together with the component 
of the director perpendicular to the layer normal, c. 

The soft mode in Eq. ([5| , denoted by Asoft can be trans- 
formed to different starting configurations of the director 
and layer normal by the following rotations 
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oft 
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^ ALftW-Q 



(8) 



where the rotation matrix Q takes the general start- 
ing layer normal ko to the z direction, and no into 
cos 6^0 z + y sin 6^0- The second rotation matrix ^ can be 

used to satisfy the requirements of the soft mode in target 
state, for example ensuring that the zx shear component 



is zero. This transformation is described in Appendix 
[a] for the case of stretching parallel to the layer normal 
in the ko = x direction, when no = cos^o^ + sin^oY- 
Although the result is analytic, the algebra is not in- 
structive, and is not presented here. The components 
of the deformation matrix for this geometry are illus- 
trated in Fig. [3J The A^^ component increases with im- 




FIG. 3. The components of the upper triangular deformation 
matrix for the Sm-C soft mode stretching parallel to the layer 
normal, with r — 2 and ^o = 0.5 radians. Initially ko = x 

and no = cos^ox + sin^oy- 

posed \xxi i-e. the sample expands in the direction per- 
pendicular to the imposed elongation. This is because 
the constraint requiring a fixed angle between the layer 
normal and director results in the director rotating into 
the z direction. The sample then expands to accom- 
modate the anisotropic chain shape, as shown in Fig. [4] 
This illustrates an unusual property of some Sm-C soft 
modes: their negative Poisson's ratio. To our knowledge 
this mechanism for negative Poisson's ratio has not been 
reported before. Alternative mechanisms of producing 
auxetic behaviour based on modifying the attachment of 
mesogens to the polymer backbone in smectic LCEs have 
been proposed and investigated experimentally [33l |34] . 
The incremental Poisson's ratio (IPR) is defined by 



d\r^ 



dXr 



(9) 



where an elongation A^^^^ is imposed and X^z is the trans- 
verse deformation. For isotropic materials, the Poisson's 
ratio must be in the range —l<u< 0.5. LCEs are 
anisotropic materials, so have Poisson's ratios outside 
this range. As the materials considered here are vol- 
ume conserving, the Poisson's ratio in the y direction is 
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FIG. 4. An illustration of the Sm-C elastomer deformation when stretching parallel to the layer normal. The director (red) 
moves out into the z direction perpendicular to the stretch axis, maintaining its tilt angle with respect to the layer normal 
(white) causing the elastomer to expand in the perpendicular direction. 
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normal, the Poisson's ratio at Xxx = 1 is given by 



When stretching parallel to the layer Eq. (12) is well founded in nematic LCEs, so serves as 
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(r — 1) cos^ 6^0 



(10) 



Substituting in typical values of ^o ^ 0.5 radians, and 
r ~ 2 for a side chain system produces v ^ —1.3. This 
is corresponds to a larger expansion than is achieved in 
auxetic foam systems [35 , albeit in only one direction. 
The extent of the soft mode in this geometry is 




■2^0, 



(11) 



hence it has no extent when ^o = 0, or when r = 1. 
Consequently the result in Eq. (10) cannot be used to 



calculate Poisson's ratio for the Sm-A phase which has 
no soft deformations. 



B. Semi-soft elasticity 

Soft modes in ideal LCEs have zero energy cost, and 
so the sample requires no force to deform. In prac- 
tise these materials have several sources of non-ideal be- 
haviour, such as compositional fluctuations and cross 
linking points that result in semi- soft behaviour. We will 
use the well known, and general form (up to quadratic 
order) of semi-soft elasticity in nematics [10] 



Fss = §a/iTr A • (I — nonj) • A^ • nn^ 



(12) 



a starting point for smectic LCEs. However, in Sm-C 
elastomers the semi-soft term in the free energy could 
in principle involve any of the directions in the problem 
including the director, and the layer normal, but we will 
neglect these effects here for consistency. 

Typical values of a are up to ~ 0.1 in nematic LCEs, 
but it may be even larger in smectic LCEs [31 . 

Some studies of semi-soft elasticity have used the fol- 
lowing simplified form 



Fss = la/j.TT 
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(13) 



which is the neo-Hookean elasticity formula. This more 
general semi-soft term gives rise to the same qualitative 
behaviour as Eq. (12). 



C. Numerical Method 



The free energy described in Eq. ([T]) , (|3| , (|4| , and ( [T2| , 

is subject to the non-linear constraints that the director 
remains of unit length and that the layer normal deforms 
as an embedded plane (Eq. [2]) . This constrained minimi- 
sation can only be performed analytically in a few cir- 
cumstances. Numerical minimisation of this free energy 
using conventional methods often results in the location 
of only local minima. We have used a simulated anneal- 
ing algorithm to minimise the total free energy, which 
finds the global free energy minimum more reliably. The 
constraint of the tilt angle of between the layer normal 



and the director can be encoded as 



n = c sin ^ + k cos 0^ 



(14) 



where the vector c is perpendicular to k. A particular 
basis is required to express c. It is convenient to use Cq, 
the starting component of no perpendicular to ko, and 
Co X ko- The vector c can be expressed as 



acos( 



b sin ( 



(15) 



where a is a unit vector constructed from the component 
of Co that is perpendicular to k, and b is perpendicu- 
lar to both a and k, i.e. b = k x a. The simulated 
annealing algorithm then minimises the free energy over 
^, and the required components of A- The global mini- 
mum derived from this was then refined using a NAG se- 
quential quadratic programming library routine. The im- 
posed constraints are implemented using Lagrange multi- 
pliers. The results of this method are in good agreement 
with the results obtained from configurations that can be 
solved analytically. 

III. ELONGATIONS OF SM-C ELASTOMERS 

We will consider four elongations to illustrate some of 
the behaviour and to build up some intuition for semi-soft 
Sm-C elastomers. The orientation of the layer normal 
and director in each case is shown in fig. [5j The mo- 
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experimental work on polydomain Sm-C elastomers [TT], 
and the mathematical studies of the Sm-C free energy 
to find which soft deformations are permitted by the for- 
mation of compatible microstructures [20^. The elastic 
behaviour of the polydomains is considerably more com- 
plicated. Understanding the elastic deformations of a 
semi-soft monodomain is a useful step towards modelling 
the experimentally more accessible polydomain sample. 
We will ignore the effect of clamping at the boundaries, 
and focus on the deformation of sheets of Sm-C elastomer 
whose mechanical properties will be dominated by the de- 
formation of the material in the middle of the long sheet. 
We will consider elongations in the x direction, together 
with the induced shear deformations. The appropriate 
deformation matrix is 



(16) 



The yx and zx components are set to zero as they would 
be resisted by counter torques. The zy component can 
be set to zero by allowing suitable rotations about the 
X axis, along which the elastomer is stretched. We will 
consider an imposed deformation A^^^^. In experiment, 
imposed stress ensembles are often used, which yield the 
same results when the stress- strain curve is monotonic. 
However, some of the stress-strain curves calculated here 
are non- monotonic, hence there are several strain values 
for a single stress value. In this case there is a difference 
between the fixed stress and fixed strain ensembles, and 
for fixed stress a Maxwell construction must be used to 
determine the strain. This is described in [36 , and briefiy 
in 3TVl 
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The model has the parameters ii^at^B^r and ^o- Typ- 
ical, ^0 ^ 30° [22], B I \i = 6 ~ 60 in well ordered sam- 
ples [T2j [151 El 7 <^^/m = c > 1 and a ~ 0.1 in smectics 
J3T| [32] , and r ~ 2 in side chain liquid crystalline poly- 
mers [1] . We will use these parameter values to illustrate 
the behaviour of the model in what follows. 
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FIG. 5. The orientation of the director and layer normal for 
each of the elongations considered, a) perpendicular to the 
layer normal (note that the director and layer normal are both 
perpendicular to the x direction initially), b) parallel to the 
layer normal, c) perpendicular to the director and d) at an 
angle '0 to the layer normal. 

tivation for these different geometries is principally the 



A. Elongation perpendicular to no and ko 

First we consider an elongation deformation in the x 
direction, with the starting layer normal ko = z and the 
starting director no = cosOqz + sin^oY, as illustrated in 
fig. [5] a). In the absence of the semi-soft term of Eq. (12) 
this deformation is as described in ^ 11 A[ The full free en- 
ergy can be minimised numerically as explained in §11 C| 
The resulting stress-strain curve, and the orientation of 
the director of this minimisation are shown in fig. [6] by 
the thick (green) lines. For the ideal Sm-C elastomer, 

this plateau ends at A^^^^ = . /-, as can be seen from the 

soft mode in Eq. ([5|. The plateau ends when the director 
has completed a rotation by 7r/2 around the layer normal. 
For non-zero values of a the onset of rotation of the layer 
normal is delayed, and it never finishes a full 7r/2 rotation. 
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FIG. 6. The stress and the angle of rotation of a Sm-C elas- 
tomer when stretched perpendicular to both the layer normal 
and the director. The model parameters are 6 = 60, r = 2, 
Oo = 0.5 (radians) and the a and c values shown on the figure. 
The thick curves (green) are from the more general numerical 
relaxation, and the black curves are calculated using the de- 
composition of the deformation matrix explained in the text. 



This is evident in the stress-strain curve, because the well 
defined stress plateau for a = becomes progressively 
less sharply defined. For a ~ 0.01 there is a pronounced 
stress plateau, but for larger values of a ~ 0.1 there is no 
plateau, merely a knee in the stress-strain curve. Fig. |6] 
also shows the effect of reducing the tilt modulus c. The 
knee in the stress strain curve becomes less pronounced, 
and the rubber hardens more slowly for larger values of 
Xxx' The retardation of the director rotation may be 
significant for piezoelectric response of these materials. 
There would be no piezoelectric response until the strain 
was above the threshold. The potential difference across 
the sample would be lower in semi-soft samples because 
the alignment of the electric dipoles associated with di- 
rector rotation is spread over a much larger deformation 
range. 

The deformation components when stretching perpen- 
dicular to k are illustrated in fig.[7| Note the sympathetic 
shears that accompany the director rotation are persis- 
tent, because the director rotation is never completed if 
a> 0. 

Numerically it is clear that with the inclusion of the 
semi-soft term there is a delay in the rotation of the di- 
rector. Some analytical progress can be made in this ge- 
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FIG. 7. The components of the deformation tensor for 
(a, 6, c, ^0, t) = (0.05, 60, 1, 0.5, 2) when stretching perpendic- 
ular to both the director and the layer normal. Note that 
the sympathetic shears persist, as the director is unable to 
complete its 7r/2 rotation. 



ometry by decomposing the deformation into three parts; 
the initial hard deformation with fixed director and layer 
spacing denoted Ahard, the soft mode Asoft and the sub- 
sequent shear and elongation after the soft mode A' |36] 



A — A • Asoft • Ahard (17) 

diag(Ai,l/Ai,l), Asoft given in Eq. (Isl), 



where Ahard 
and 



C -q 

A' = I 1/c 

1 



(18) 



This deformation matrix can be substituted into the free 
energy terms of Eq. (IT]), (Is]), and ([12]) (assuming that 
c ^ oo, so that d = OoY The problem is then reduced to 
a minimisation over the variables XiX^V and ^, with the 
constraint that the total Xxx is prescribed. The threshold 
before the onset of director rotation can be calculated by 
setting C = 1 and f] = 0^ then performing a series expan- 
sion of the free energy in soft mode rotation angle (f). The 
leading term is (9(^^), and when this term becomes neg- 
ative a non-zero value of (j) will lower the free energy. To 
leading order in (Ai — 1), this coefficient becomes negative 
when Ai is approximately 

Ai= 1 + Sr^a/{1 y 29r - 29r^ - r^ y ra y 3br^a 
+4r^acos2l9 + (r - l)((r - 1)^ + ra) cos 4(9) (19) 

This value is slightly smaller than the corresponding 
threshold to director rotation in nematic elastomers of 
-^1 = r-i-ar i'^^' Intuitively this is because in the Sm- 
C phase the deformation is restricted to two dimensions 
by the layer spacing constraint. Consequently there is a 
larger contraction in the direction perpendicular to the 
stretch which causes the elastic free energy to rise faster, 
and hence the director rotation to start earlier in Sm-C 
LCEs as compared to the nematic phase. 



The minimisation of the free energy over XiX^V ^^^ 
(f) produces results that are in good agreement with the 
more general numerical method. The results are shown 
by the black lines in Fig. [6] 



B. Elongation parallel to kg 

Elongation parallel to the layer normal is illustrated in 
Fig.lslb). The initial layer normal is given by kg = x and 
the director no = x cos ^o + z sin ^o- Using the form of de- 
formation matrix described in Eq. (16), the free energy 



can again be minimised using the numerical technique 
described in 811 CI The results for various values of the 



semi-soft parameter a are illustrated in fig.[8] For c ^ oo 
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FIG. 8. a) The stress-strain response for a semi-soft Sm-C 
elastomer stretched parallel to the layer normal. The model 
parameters are 6 = 60, r — 2 and ^o = 0.5, and the values of 
(a, c) shown in the figure, b) The corresponding IPRs for the 
stress-strain curves. 

the first part of the stress-strain curve is determined by 
the smectic layer modulus B. The semi- soft term pre- 
vents the rotation of the director, and the layer spacing 
increases. Once the force required to increase the layer 
spacing is comparable to that required to rotate the di- 
rector the semi-soft mode begins. The stress-strain curve 
has negative slope once director rotation starts. As ex- 
plained in §11 A| there is a negative IPR in this geometry 
as the director rotates around the layer normal into the 
direction perpendicular to the stretch axis (see Fig.jsjb)). 
This lateral expansion, combined free energy expression 
for the semi-soft elasticity, results in the negative stiff- 
ness. For larger values of a the Poisson's ratio becomes 



less negative. 

The rotation of the layer normal and director, and the 
deformation components are illustrated in Fig.|9] The ex- 
pansion of the sample in the z direction is clearly visible 
at the onset of rotation, as are the usual shear compo- 
nents that accompany a soft mode. For finite values of 
c the deformation becomes more complicated; before the 
threshold the director rotates towards the layer normal 
and the sample shears, which itself results in movement 
of the layer normal. There is both an increase in the 
threshold to the start of rotation, and a reduction in the 
amplitude of the semi-soft deformation. This is because 
the shearing before director rotation results in rotation 
of the layer normal, and there is a reduction in the tilt 
angle before the onset of shearing. 




1.16 1.18 



•^xx 



FIG. 9. For the parameter values (a, c) = (0.05, oo) and 
6 = 60, r = 2 and ^o = 0.5 radians for stretching parallel 
to the layer normal, a) Shows the director and layer normal 
rotation, b) the shear components, and c) the diagonal com- 
ponents of the deformation tensor when stretching parallel to 
the layer normal. 



The soft mode in this geometry can be calculated ana- 
lytically, as explained in ^ 11 A[ and its amplitude is given 
in Eq. (11). Some analytic results can be obtained by 



decomposing the deformation as follows. 

A = A • Asoft • Ahard • 

Both A' and Ahard have the form 

r, C 
C 
l/(r?0 



(20) 



(21) 



The appropriate soft mode must be calculated based on 
the rotated layer normal, as the shear component ( will 
cause it to rotate. 

The onset of director rotation can be obtained by sub- 
stituting back into the free energy, expanding in terms 
of (p up to quadratic order. When the coefficient of the 
O(0^) term is negative, the soft mode becomes active. 
This happens when 



Ai 



2 



3ar^ 



b{r 



1 ) ^ cos^ ^0 4p^ 6^ c cos^ ^0 



-\a') 



(22) 



Note for smaller values of c this is inaccurate because the 
shear is only expanded up to quadratic order. The oc- 
currence of a and b in this expression correspond to the 
competition between the stretching of the layer spacing, 
and the semi-soft elastic term keeping the director fixed 
in the matrix. The threshold predicted by this calcula- 
tion is consistent with the numerical results for large b 
and c. 



1. Scalar model of negative slope region 

The unusual response above for the Sm-C soft mode 
can be illustrated for a much simpler deformation. Con- 
sider an elongation with a diagonal deformation matrix 
of an imposed \xxi ^zz given by 



AiX 



A 
4' 



(23) 



with \yy determined by volume conservation. The pa- 
rameter A here controls the initial rate of expansion of 
the material. Its Poisson's ratios are —A, and l + A. This 
is similar to the soft mode in a Sm-C illustrated in Fig.|3] 
The deformation in Eq. (23) can be substituted into a 



neo-hookean model such as Eq. (13), which is broadly 



similar to the semi-soft elastic energy term. The result- 



ing stress-strain curve is shown in fig. 10 It can be seen 
from this plot that for sufficiently large values of A the 
stress-strain curve has a negative slope similar to stretch- 
ing the Sm-C LCE parallel to the layer normal. For some 
geometries the Poisson's ratio is sufficiently negative to 
result in a negative stiffness. The configurational entropy 
of the perpendicular degrees of freedom decreases as the 
sample expands resulting in a positive contribution to the 
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FIG. 10. For the scalar model of the negative stress strain 
curve described in the text a) shows the stress-strain curves 
for A — 0.5,0.75,1, and b) the deformation components for 
A—\ for the scalar auxetic model. 



stress. Once lateral expansion starts to slow sufficiently 
there is a weaker contribution to stiffness of the sample 
from the perpendicular degrees of freedom and the stress 
starts to drop, which produces a negative slope in the 
stress-strain response. By tuning the parameter A in the 
model, the balance between the parallel and perpendic- 
ular degrees of freedom can be altered, and the stiffness 
changed from negative to positive. 

This scalar model shows that the negative stiffness is 
a result of the lateral expansion during the Sm-C soft 
mode, and not due to the form of the semi-soft elastic 
term. 



C. Elongation perpendicular to no 

Stretching perpendicular to the initial director, no is 
illustrated in fig. [sjc). The results for the numerical cal- 
culation of the stress-strain curve for this geometry are 
shown in Fig. [TT] This geometry has the remarkable fea- 
ture that Vzz -^ —CO when a — > 0, as shown in Fig. 11 
For larger values of ol the Poisson's ratio becomes less 
negative. The jump in the director also causes a discon- 
tinuity in the IPR, and a sudden increase in the width 
of the sample. Note that in this geometry there is a 
discontinuity in the stress-strain curve, in addition to 
the negative stiffness. The discontinuity in the stress- 
strain curve is accompanied by a jump in the director 
as shown in Fig. [12] Intuitively the discontinuity arises 
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FIG. 11. a) The stress-strain curves for stretching perpendic- 
ular to the layer normal, b) the IPR for various parameters 
6 = 60, r = 2 and ^0 = 0.5 and values of (a, c) shown on the 
figure. 



because when the director jumps the long axis of the 
polymer shape tensor jumps towards the elongation di- 
rection. Consequently the natural length of the rubber 
in this direction is increased, so there is corresponding 
drop in the stress. 

The jump in the director can be understood from the 
properties of the soft mode in this geometry. We can 
approximate the first part of the total deformation (until 
the end of director rotation) as a hard deformation where 
there is no director rotation, followed by a soft mode 



A — Asoft • Ahard • 



(24) 



The soft mode in this geometry can be calculated ana- 
lytically as explained in appendix [Aj Whilst its analytic 
form is algebraically very long, the amplitude of the soft 
mode has a much simpler expression, and is given by 

Kx= (3 + r(7r - 2) + 4(r^ - 1) cos 2(9o+ 

(1 + (2 - 3r)r) cos4l9o)^/^ /(2V2p). (25) 

The hard part of the deformation has only diagonal ele- 
ments, and an xz shear component. 



A 



Ahard — 



^xx 


K, 





1/(A,,A,,) 





K. 



(26) 



Substituting this into the full free energy density yields 
an approximate solution to the minimisation problem. 
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FIG. 12. When stretching perpendicular to the director, for 
the parameter values (a, c) = (0.05, 00) and 6 = 60, r = 2 
and ^0 = 0.5 radians a) shows the director and layer nor- 
mal rotation, b) the shear components, and c) the diagonal 
components of the deformation tensor when stretching per- 
pendicular to no. 



where the director rotation is assumed to be continuous. 
The free energy density in this case is shown in Fig. 13 
The analytic solution with continuous director rotation 
has higher free energy for the first part of the deforma- 
tion. Hence, the elastomer initially stretches without di- 
rector rotation. If the director were to start rotating, 
then the form of the soft mode results in rapid rotation 
of the director, and an infinite slope in the free energy. 
However, the rate of increase slows, and eventually the 
state with a rotated director is lower in free energy than 
that with a fixed director. At this point the director 
jumps to the new orientation. There is a discontinuity in 
the slope of the free energy at this point, or equivalent ly 
a jump in the stress. 

This behaviour is not solely a result of the semi-soft 
energy term, but again is a result of the shape of the soft 
mode, combined with a general semi-soft elasticity term. 
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FIG. 13. The free energy calculated numerically, and the free 
energy trajectory of the semi-soft mode with continuous di- 
rector rotation when stretching perpendicular to the director. 
Here b — 60, r — 2^0q — 0.5 and (a, c) are shown on the figure. 



These calculations are based on an equilibrium model 
of a Sm-C elastomer. In practise kinetic terms, such as 
viscosity would smooth out the sharp jump demonstrated 
here. 



1. Scalar model describing stress discontinuity 

The semi-soft behaviour of Sm-C elastomers is char- 
acterised by two deformation modes; before the onset of 
director rotation, and afterwards. A scalar model that 
exhibits the same behaviour when stretching perpendic- 
ular to the director can be developed based on represent- 
ing each of these deformation modes as a spring, and de- 
forming the two springs in series. The total strain is the 
sum of two deformation modes corresponding to keeping 
a fixed director eu, and rotating the director esM 



ex = eu + esM- 



(27) 



The two modes of deformation have different energy 
penalties, the first arises from a simple uniaxial defor- 
mation, so in a neo-hookean energy model will result in 
a free energy term of the form 



1 ^ ^2 



(28) 



where Ki corresponds to the shear modulus of the rub- 
ber. The second arises from the soft mode, which has 
a singular edge in the contraction of the rubber as it is 
stretched. The zz component in the soft mode is ini- 
tially of the form \zz = 1/(1 + {^xx — 1)^) (where here 
^xx — 1 = esM- When this is put into the neo-hookean 
free energy, it results in free energy terms to leading order 
in esM of the form 



^SM 



2^2esM 



(29) 



where K2 is the corresponding shear modulus for this 
mode. In the case of the semi-soft Sm-C elastomer, this 



term arises because of the rapid rotation of the director 
during the start of the soft mode. 
The total free energy is then 



Ft = iif (ex - esM)' + h^: 



a^SM' 



(30) 



where first spring in this system is hookean, and the sec- 
ond is non-linear, being infinitely stiff at zero strain for 
< /3 < 1, but softening rapidly as strain increases. 
This should be minimised over esM to determine the dis- 
tribution of strain between the two springs. It can be 
solved analytically for /3 = 0.5. The behaviour of this 
model is illustrated in Fig. 14 For small (3 this system 



has a discontinuity in the stress-strain curve, but as /3 
is increased the stress-strain response becomes continu- 
ous. The free energy as a function of esM is also illus- 
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FIG. 14. a) An illustration of a discontinuous stress-strain 
curve for the scalar model described in the text, b) the free 
energy as a function of the variable esM for fixed total strain 
values. Here Ki — 10 and K2 = I. 



trated in Fig. 14 For small values of ct there is only 
one minimum at esM = 0, corresponding to no strain of 
the second spring. However, as the total strain increases, 
the second mode of deformation becomes activated and 
there is a minimum for larger values of esM- Since there 
is a barrier between the two minima, the transition is 
first order, so there is a jump in the equilibrium value of 
esM- For larger values of /3 the phase transition becomes 
continuous, and the stress- strain curve no longer exhibits 
a jump. 

This behaviour is analogous to that of the semi-soft 
Sm-C elastomer as the free energy exhibits a discontinu- 
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ity when stretched perpendicular to the director (where 
the soft mode has a singular edge). Larger values of /3 
correspond to stretching at a larger angle to the director, 
where the soft mode does not have such a rapid rotation 
of the director, and a corresponding sharp drop in the 
lateral dimension. If the angle between the director and 
the elongation direction is large enough, then the stress- 
strain response becomes continuous as we will see in the 
next section. 



D. Elongation at an angle ip to the layer normal 

The last deformation we consider is shown in fig. pld). 
The numerical solution of stress-strain curve associated 



a) 



with this geometry is shown in Fig. [15] The stress-strain 
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FIG. 16. For stretching at an angle of -0 = 0.65 radians to the 
layer normal a) shows the director and layer normal compo- 
nents in the xy plane, b) the shear components of the deforma- 
tion, and c) the diagonal components of the deformation for 
stretching at an angle of -0 = 0.65 radians to the layer normal, 
for the case (a, c) = (0.05, 00), and b = 60, r — 2^0q — 0.5. 



FIG. 15. a) The stress- strain curves for stretching at an angle 
of '0 = 0.65 radians to the layer normal, for 6 = 60, r = 2, ^0 = 
0.5 radians and various parameter values (a,c), and b) the 
Poisson's ratio in this geometry. 

curve is continuous in this geometry, but again has a 
pronounced negative slope. There is a negative IPR of 
^ — 1.5 that is roughly independent of the semi-soft pa- 
rameter. The expansion of the sample that accompanies 
the rotation of the director can be seen in Fig. [16] 



IV. DISCUSSION 

The first three deformations considered above in Fig.jSJ 
(a-c), when made with clamped boundary conditions, 
would not be soft even without the semi-soft elastic term. 



This is because no microstructure can be constructed 
from the soft deformations that is compatible with the 
boundary conditions, due to the shear components in 
the Sm-C soft mode [20 . However, the properties of a 
long sheet of Sm-C LCE may approximate this behaviour 
as the centre of the sample could deform without rigid 
boundary conditions. The final deformation in Fig. Isld) 
can be performed with clamped boundary conditions in 
the soft case. In the semi-soft case the sample starts to 
shear before the onset of rotation, which is not compati- 
ble with clamped boundaries, so in experiment it may be 
even stiffer initially due to this additional constraint on 
its deformation. 

The maximum lateral expansion can be deduced from 



the soft mode presented in ^11 A The shear components 
are transformed, through a rotation, into an elongation. 
At ^ = 7r/2 the maximum lateral expansion occurs (in 
the y direction for the example given in the text), and 
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has a value of 

The region of negative slope in the constitutive models 
reported here is typically explained by a Maxwell Con- 
struction. Similar behaviour occurs in the Van der Waals 
gas model which has a region of negative slope in the 
pressure-volume curve. Here there is a two phase region 
consisting of a mixture of the liquid and gas phases. In 
solids the two deformations on either side of the insta- 
bility must be compatible to form a mixture [37] • The 
system should then disproportionate, adopting a mixture 
of the two deformations to achieve the externally imposed 
strain. The first order type phase transition seen in the 
example stretching perpendicular to the layer normal can 
result in hysteretic behaviour as the system jumps from 
one energy well to another. The rate of the deforma- 
tion in comparison to the sample relaxation times may 
also result in hysteresis [38 . There is interest in negative 
stiffness materials [39] for applications such as sealants, 
stiffening composites, and creating meta-materials hav- 
ing a negative refractive index to sound waves. 

Experimental work reporting mechanical testing on 
Sm-C monodomains has not been reported. Whilst it is 
anticipated that these monodomains should exhibit soft 
elasticity, the addition of the semi-soft elasticity term to 
the model suggests that these effects may be difficult to 
observe for large semi-soft parameter a. When stretching 
perpendicular to both the layer normal and the director, 
the semi-soft term may prevent any stress plateau being 
observed, instead only a shoulder is visible in the stress- 
strain response. 

Although we have only considered the deformations of 
monodomains here, the results inform model predictions 
for poly domains. Poly domains are difficult to model be- 
cause of the requirement of ensuring that adjacent do- 
mains deform in a compatible way. A simplifying ap- 
proximation used to model a polydomain is to assume 
that it consists of an array of monodomains that deform 
at the imposed external strain, but are independent from 
each other. If we deform the pseudo-monodomain shown 
in Fig. [T] by stretching in the x direction, then defor- 
mation component \yy averaged over all the domains is 



illustrated in fig. 17 for 50 domains. This figure shows 
that there is a negative IPR as the director in each of the 
domains jumps causing them to expand. The curve illus- 
trated here is jagged because the alignment of each do- 
main jumps at a slightly different threshold. The expan- 
sion of the film thickness, and the energy loss as a result of 
the jump in the director orientation in this geometry may 
be observable in experiments on pseudo-monodomains 
pT| |4Q . The larger values of deformation reported in ex- 
periment before the knee in the stress-strain curve point 
to a much larger value of a than in the illustrative plot 



in Fig. 17 




The features of the smectic-C model described here 
would be present in a wide range of models that have 
soft modes of nematic elastomers but incorporate the 
constraint on the director to remain at a fixed angle 
to the layer normal. However, validation of these mod- 
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FIG. 17. The average value of \yy for 50 domains in a pseudo- 
monodomain illustrated in fig. [l] as a function of \xx assum- 
ing they all experience the same strain, and deform indepen- 
dently. Model parameters are b — 60, r = 2,^o = 0.5, c = 
oo, a = 0.05. 



els await either experimental work on mechanical testing 
of Sm-C monodomains, or theoretical work on pseudo- 
monodomains to link up with existing mechanical exper- 
iments on pseudo-mododomains. 



V. CONCLUSION 

We have studied a model of monodomain Sm-C LCEs 
with the inclusion of a semi-soft elastic term to describe 
imperfections in the elastomer. As result of the negative 
incremental Poisson's ratio inherent in the soft modes 
of a Sm-C monodomain, the mechanical properties of 
a semi-soft monodomain are unusual. When stretching 
perpendicular to the layer normal and the director, the 
response is reminiscent of a nematic elastomer. A finite 
force is required to deform the LCE and initial the ro- 
tation of the director. However, the stress plateau is 
less well defined for larger values of semi-soft param- 
eter a; it is reduced to a shoulder in the stress-strain 
response. When stretching parallel to the layer normal 
the elastomer again exhibits a threshold to director ro- 
tation. Once director rotation has started the elastomer 
has a negative incremental Poisson's ratio, and a negative 
stiffness. A negative incremental Poisson's ratio of up to 
V r^ —1.5 has been found for typical model parameters. 
This arises because the director rotates in a direction per- 
pendicular to the stretch axis due to the constraint of the 
layer normal. This more detailed understanding of mon- 
odomain deformations of Sm-C elastomers might prove 
useful in understanding recent mechanical and piezoelec- 
tric experiments on polydomain Sm-C elastomers. 
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Appendix A: Transforming the soft mode for 
different starting configurations 

The soft mode given by Eq. ([5| can be transformed to 
other geometries by a pair of rotation matrices. For ex- 
ample, consider the case of a Sm-C elastomer stretched 
parallel to the layer normal. Let us assume that start- 
ing layer normal is ko = x and the starting director is 
no = cos^ox + sin^oy- The soft mode for this configu- 
ration that is an upper triangular matrix, as described 
in Eq. (16), can be found as follows. From the reference 



configuration a body rotation is performed such that the 
layer normal ko is parallel to the z axis. In this case, a 
90° rotation about the y axis 



Q = 




(Al) 



After this rotation the director is given by n = z cos ^o + 
ysin^o- Note that in general an additional rotation 
around the z axis may be required to ensure the direc- 
tor is in this orientation. This is the initial configuration 
for the soft mode given in Eq. ([s]). The director now ro- 



tates by an angle <j) around the new layer normal, and 
the sample executes the soft mode. Finally a rotation of 
the target state is performed such that the deformation 
matrix has the form described in Eq. (16). This rotation 



matrix is in general simpler if we first undo the rotation 
Q. The rotation matrix P_ is described by three angles: 



cos \l)z sin \Ijz 

— sin ipz cos ipz 

1 

cos iljy sin iljy 

1 

— sin ipy cos ipy 



1 

cos ipx sin ipx 
— sin ipx cos ipx 



(A2) 



The three angles V'x, ^y and tp^ can be calculated by sub- 
stituting into the equation 



\ = P-' 



r-x. 



Asoft • ' 



(A3) 



and ensuring that the three lower triangular elements of 
A are zero. In the case of stretching parallel to ko and 

perpendicular to no the soft mode can be calculated an- 
alytically. The algebraic expressions for these soft modes 
is long, and unedifying so will not be presented here. 
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